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ABSTRACT 

We have used the two-degree field (2dF) multi-fibre spectrograph of the Anglo- 
Australian Telescope to search for candidate members of the unusual globular cluster 
uj Centauri at and beyond the cluster tidal radius. Velocities with an accuracy of ~10 
km s _1 were obtained for 4105 stars selected to lie in the vicinity of the lower giant 
branch in the cluster color-magnitude diagram and which cover an area on the sky 
of ~*2.4 x 3.9 deg 2 centered on the cluster. Within the velocity interval 190-270 km 
s -1 , the cluster member candidates have a steeply declining surface density distribution 
consistent with the adopted cluster tidal radius of 57'. For the stars in the sample 
beyond the tidal radius, an analysis of line-strengths from the spectra, as well as radial 
velocities, identifies only six stars as possible candidates for extra-tidal association with 
the cluster. If all six of these stars are indeed related to the cluster, then a maximum 
of 0.7 ± 0.2 per cent of the total cluster mass is contained in the region between one 
and two tidal radii. Given this limit, we conclude that there is no compelling evidence 
for any significant extra-tidal population in u Cen. The effects of tidal shocks on the 
outer parts of the cluster are consistent with this limit. Theories for the origin of lo Cen 
frequently suggest that the cluster is the former nucleus of a tidally stripped dwarf 
galaxy. Our results require that the stripping process must have been largely complete 
at early epochs, consistent with current dynamical models of the process. The stripped 
former dwarf galaxy stars are therefore now widely distributed around the Galaxy 
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Introduction 



The stellar system uj Centauri is unusual in many respects. When considered as part of the 
Galaxy's globular cluster population it is among the most luminous and the most massive, with 



estimates for the cluster mass ranging from 2 to 5 x 10 6 Mq (e.g. van de Ven et al. 



2006 



Meylan 



et al. 1995). Unlike other globular clusters though, uj Cen exhibits the characteristics of system 



that has undergone substantial chemical self-enrichment, with its member stars showing a wide 
spread in abundance (e.g. Freeman Rodgers||1975 Norris Da Costa||1995 Pancino et al.|[2 002), 



and complex distributions of element-to-iron abundance ratios (see Romano et al. 2007 and the 
references therein). The abundance picture is further complicated by the discovery that the lower 



main sequence of the cluster is bimodal (Bedin et al. 2004), with the bluer sequence having a 



higher abundance (Piotto et al. 2005). These observations can be explained if the metal-rich but 
bluer sequence has a substantially higher helium abundance (Y ~ 0.40, see: Norris|[2004 Piotto et 



al. 2005), but generating such large helium abundances within the context of chemical evolution 



models is difficult (e.g. Romano et al.|20"07 ). An extensive spectroscopic survey of uj Cen post-main- 
sequence members that characterizes the many and diverse properties of this cluster is presented 



m 



van Loon et al. (2007) 



The dynamics of uj Cen are also unusual. Unlike most globular clusters, it has a notably 



flattened shape and rotates comparatively rapidly (see van de Ven et al. 2006 and the references 
therein). Moreover, the stellar abundances and stellar kinematics are correlated, with the more 



metal-rich stars being not only more centrally concentrated but also kinematically cooler (e.g. Norris 
et al.|199~ Sollima et al.|2005 ), although see also Pancino et al. (2007). Possible substructures with 



the cluster may also exist. For example, van de Ven et al. (2006) find evidence in their kinematic 



analysis for a separate disk-like component between 1 and 3 arcmin from the cluster center. Finally, 
the orbit of uj Cen is also unusual in that it is tightly bound and retrograde with apo- and peri- 



galactocentric distances of 6.2 and 1.2 kpc, respectively, and a period of 120 Myr (Dinescu et al. 



1999b). The cluster never rises very far above the Galactic plane: Dinescu et al. (1999b) give z. 



as 1.0 ± 0.4 kpc. 



These characteristics have led to increasing support for the idea (cf. Freeman 1993 Norris et 



al. 1996 1997) that uj Cen has not evolved in isolation, but in fact represents the nuclear remnant 



of a now disrupted nucleated dwarf galaxy that has been accreted by the Milky Way. Models of the 
orbital decay and disruption of a nucleated dwarf, which is presumed to lie initially far from the 
Galaxy's centre, have shown that this process is dynamically plausible, with most of the dwarf's 



mass lost as the nucleus settles into an orbit like that of the present-day uj Cen ( Bekki & Freeman 



2003 Tsuchiya et al.|2003 ). These dynamical models, however, do not cover a full Hubble time and 
while the likely present-day outcome is that the surviving nucleus is an isolated bound remnant, 
one may seek to confirm this prediction observationally. For example, in the model of Bekki & 



Freeman (2003), the mass of the stellar envelope surrounding the nucleus, which comes from the 
disrupting dwarf, is still fully ~50% of the mass of the nucleus at the end of the calculations 2.6 Gyr 
after the initiation of the interaction. At this point the nucleus has achieved an orbit comparable 
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to that of oj Cen. The most likely signature of any remnant from such an interaction would be in 
the form of structure beyond the tidal radius of the cluster (nucleus) , probably in the form of tidal 
tails. 

Extra-tidal structure, particularly in the form of tidal tails, is well established around a number 
of globular clusters (e.g. Grillmair et al.fl995 Leon et aT||2000 Grillmair &z Johnson||2006 and the 



references therein) , with the most striking example being Pal 5 ( Odenkirchen et al.||2003 Grillmair 
Dionatos||2006[ ). The origin of these features though lies with the interaction between the cluster 
and the tidal field of the Galaxy rather than as a remnant from any merger process. Tidal tails 



have been claimed to exist in oj Cen by Leon et al. ( 2000 1 who analyzed the spatial distribution of 



color-magnitude diagram selected stellar samples, chosen to maximize the contrast between field 
and cluster stars in their photographic data. For oj Cen they found "two large and significant tidal 
tails" extending up to ~2 degrees to the north and south of the cluster and estimated that the 
features contained up to 1% of the cluster mass. They noted, however, that dust absorption, which 



they did not account for, might affect the star counts. Subsequent work by Law et al. (2003) (see 



also Coleman & Da Costa 2008), however, showed that these features are much less prominent 



when reddening variations are properly accounted for. 

In this paper we present the results of an extensive spectroscopic survey in which we use radial 
velocity and line strength information to investigate the number and distribution of oj Cen stars 
both within and beyond the cluster tidal radius, and to place constraints on the outer structure of 
the cluster and on any extra-tidal population present. The following section describes the selection 
of the sample, the observations and the derivation of the radial velocities. Section 3 provides the 
analysis for stars both in the cluster and outside the cluster tidal radius. The main results are 
discussed in §4 and the conclusions are summarized in §5. 



Observations and Reductions 



2.1. Sample Selection 



In a forthcoming paper, Coleman Sz Da Costa (2008) will present the results of an imaging 
survey of a region ~2.4 x 3.9 deg 2 in size centered on oj Cen. The V and /-band photometry 
was obtained with the 8k x 8k CCD mosaic camera mounted on the ANU lm telescope at Siding 
Spring Observatory. The observation, reduction and analysis techniques employed are essentially 



the same as those in Coleman et al. (2005a) and Coleman, Da Costa Sz Bland-Hawthorn (2005b). 



Fig.[T]shows a color-magnitude diagram (CMD) from the survey for stars that lie in the radial range 
10' < r < 25' from the cluster center. The principal cluster sequences such as the horizontal and 
giant branches, and the main sequence turnoff, are readily visible although equally, considerable 
field contamination is visible, a consequence of oj Cen's relatively low galactic latitude (b ~ +15°). 
It is fortunate, however, that the radial velocity of the cluster is relatively large (+232 km s _1 ) 
permitting a generally clean separation of members from non-members via the measurement of 
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radial velocities (e.g. Stanford et al. 20061. 



A substantial radial velocity survey is required, nevertheless, in order to generate a sample of 
likely members in the outer parts of the cluster and beyond. The two-degree field (2dF) multi- 



fibre instrument on the Anglo- Australian Telescope (Lewis et al. 2002) provides an efficient way 
to conduct such a survey. The 2dF is a double-buffered system allowing the fibres for the next 
field to be configured while observing the current field. The reconfiguration time is of order 45-60 
minutes. Consequently, the most efficient survey rate is achieved by selecting a magnitude interval 
that allows spectra with appropriate signal-to-noise for velocity and line-strength measurements 
to be obtained in total exposure times comparable to the reconfiguration time, for a given 2dF 
instrumental setup. The instrumental setup chosen used the 1200B gratings centered at A4250A 
giving a wavelength coverage of approximately AA3700-4800A at a scale of 1.1 A/pix. The resulting 
spectra have a resolution of ~2.5A. These considerations led to the choice of the magnitude range 
15.4 < V < 16.75 for the initial sample selection. The corresponding color limits were then set 
by allowing for photometric errors and by noting that the giant branch of u Cen has a significant 
intrinsic color width. The selection box is outlined on Fig. [TJ We note that over the area of the 



imaging survey there are significant variations in the amount of interstellar reddening (cf. Leon et 



al. 2000 Coleman & Da Costa 2008). Consequently, the limits of the selection box were altered 
by the appropriate amounts given the difference in mean reddening between each CCD mosaic 
field and the central field. The largest shift applied was 0.08 mag in V — I and 0.18 in V . The 
initial selected sample was then modified to exclude any candidate that lay within ~10" of a bright 
neighbor. 

The imaging survey covers a rectangular region ~2.4° (EW) by 3.9° (NS) in size while the 2dF 
field is circular. We therefore chose four 2dF field centers to the N, S, E and W of the cluster center 
to maximize the coverage of the imaging survey area. Stars closer to the cluster center than 20' 
were excluded to minimize crowding concerns while still leaving enough likely members to define 
radial velocity and line strength distributions. Six configurations were generated for both the N and 
S field centers and four for the E and W fields with a total of approximately 5500 candidates. For 
each field approximately 15 stars were kept in common between all configurations while a modest 
number of stars in the overlap regions between fields were retained in the configurations for both 
field centers, all to allow monitoring of any systematic effects. Each configuration also included 
~40 fibres allocated to blank sky regions, evenly distributed between the two 2dF spectrographs. 
These regions were selected from the imaging survey to contain no detected stars (i.e. V > 20) 
within a radius of 10". 



2.2. Observations 

The observing run took place over 5 nights in 2003 May and a total of 16 configurations were 
observed: six in each of the N and S fields and two each in the E and W fields. The observing 
sequence for each configuration consisted of three 1800 sec exposures preceded and followed by arc 
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Fig. 1. — A color-magnitude for stars with radial distances from the center of oj Cen in the range 
10' to 25' taken from the photometry of Coleman & Da Costa (20081. The region outlined on the 
cluster giant branch defines the selection box for the candidate sample observed spectroscopically. 
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lamp exposures and a fibre flat field exposure. Multiple exposures, including exposures recorded in 
each of the two 2dF spectrographs, were also obtained of six bright radial velocity standards. These 
spectra provide templates for cross-correlations with the program star spectra and for setting the 
velocity zero point. 

The observed data sets, which come in pairs as the 400 2dF fibres feed two separate spectro- 
graphs each imaging 200 fibres onto a CCD detector, were reduced with the 2dF reduction program 
2dfdr. The fibre flat field exposure was used to generate the 'tram line' map, which was then ro- 
tated and shifted to track the pixel locations of the individual fibre spectra. The fibre spectra were 
then 'fit' extracted with the wavelength calibration coming from the appropriate arc lamp expo- 
sures. The fibre spectra for the blank sky regions were then median-combined to form a single sky 
spectrum, which was then subtracted from each object spectrum using a weighting appropriate for 
the transmission of the object fibre. The relative fibre transmissions were determined either from 
twilight sky exposures with the same configuration, or from median combined offset sky exposures 
obtained as part of the observing sequence. The reduced spectra from the individual exposures for 
each configuration were then median combined to remove cosmic-ray contamination. Examples of 
the reduced spectra are shown in Fig. [2j 



2.3. Radial Velocities 

The radial velocities were determined by cross-correlation techniques using the IRAF routine 
fxcor. The velocity analysis was carried out independently for each of the two 2dF spectrograph 
data sets, including the use of template spectra observed with the appropriate spectrograph. Only 
spectra that had at least 400 counts at the G-band (~A4300A) were included in order to minimize 
the velocity errors. The primary template spectra employed were those of the V=6.20 F8 IV-V 
star HD162396 which provided the best spectral match to the program stars. The correlation was 
carried out over the wavelength range AA3800-4700A, a spectral region that encompasses a number 
of strong features such as the Call H and K lines, the G-band of CH and the H<5 and H7 hydrogen 
lines (cf. Fig. [2). 

Heliocentric corrections were then applied to correct for the Earth's motion resulting in veloci- 
ties for 4105 program objects relative to the velocities for the two template spectra (i.e. the velocity 
zero points for the two data sets) which remain to be determined. These were determined through 
the observations of standard stars of known velocity: the zero points are set so that the mean 
difference between the observed and catalogue velocities for the standards is zero. Six observations 
of four standards in addition to that of the primary template were employed for the spectrograph 1 
data set while five observations of the same four standards plus that of the primary template were 
used for the spectrograph 2 data. The uncertainty in the zero points is then the standard error of 
these mean differences. The values are ±2 km s~ x for the spectrograph 1 data and ±3 km s^ 1 for 
the spectrograph 2 data. The overall process was verified by making use of the eight programs stars 
that had at least two velocity determinations in both the spectrograph 1 and the spectrograph 2 
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Fig. 2. — Examples of the reduced spectra. The uppermost spectrum is typical of the low velocity 
relatively metal-rich field star contaminants. The middle two spectra have radial velocities com- 
patible with uj Cen membership and the stars lie well inside the cluster tidal radius. The lowermost 
spectrum is for a star with a radial velocity compatible with that for the cluster but which lies well 
outside the tidal radius. 
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data. The mean velocity difference (spectrograph 1 - spectrograph 2) for these stars is 2 km s _1 
with a standard error of the mean of 4 km s _1 , indicating excellent consistency between the two 
data sets. The errors in the individual program star velocities were then determined from the 19 
program stars that had at least five observations during the observing run, including observations 
across both data sets. The mean standard deviation from these sets of multiple observations, which 
we take as the error associated with a single program star velocity determination, is 11 km s _1 ; 
the uncertainty in this value is ±1 km s -1 . Fig. [3] shows the spatial location of the 4105 program 
stars with measured velocities. 



3. Analysis 



3.1. Stars within the tidal radius 



Figure |4] shows velocity histograms for all stars inside the adopted cluster tidal radius of 57'. 



This value has been taken from the compilation of Harris (1996). Given the importance of this 



parameter, a brief discussion is appropriate. Trager et al. (1995) list a core radius of 2.6' for the 



cluster and a concentration index c=1.24, which implies a tidal radius of 45'. Da Costa (1979) 



derived similar values (2.6', 1.22, 44') from King (1966) model fits, but noted that the observed 



data indicated a larger value for the tidal radius. The best-fit dynamical models of Meylan (1987) 



and Meylan et al. (1995), which are a sizeable improvement over the simple King (1966) models 



have tidal radii between 46' and 60', with the larger values showing a better fit to the outermost 
observed surface density data points. The data presented in §3.3 also imply a tidal radius consistent 



with the Harris ( 1996 ) value. Fortunately, none of the results of this paper depend significantly 



on the actual value of the cluster tidal radius, provided it is similar that adopted to within a few 



arcmm. 



The histogram in Fig. [4] outlined by the solid line shows the velocity distribution for the 
stars with radial distances from the cluster center between 20' and 30'. It is clearly bimodal with 
one population, predominantly foreground disk stars, centered near zero velocity and a second 
population, predominantly cluster members, centered near the cluster velocity of 232 km s _1 (e.g. 



van de Ven et al. 2006). The mean velocity for the 119 stars within the velocity interval 190-270 



km s _1 (chosen to encompass all likely members given the distribution in Fig. |4| is 229 km s _1 
with a dispersion of 14 km s . The mean velocity agrees well with the known cluster velocity 
given the uncertainty in the velocity system zero point discussed above. Moreover, the observed 
dispersion value is consistent with that expected from the combination of the intrinsic dispersion 
of the cluster at these radii, ~8 km s _1 , (cf. van de Ven et al. 2006) and the individual velocity 



errors. These stars will be referred to as the inner cluster sample. 

The histogram outlined by the dot-dash line in Fig. [4] is for stars with radial distances between 
30' and 57'. Because of the larger area and the decreasing density of cluster stars with increasing 
radial distance, the contaminating population now considerably outweighs the cluster population, 
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Fig. 3. — The location on the sky of the 4105 program stars with measured radial velocities. The 
rectangular region outlined with dashed lines is the extent of the imaging survey described in 



Coleman & Da Costa (20081 while the four circles marked with dotted lines outline the four 2dF 



fields. The center of the cluster is labeled with a circled plus-sign and the red circle marks the 
adopted cluster tidal radius of 57'. No stars were observed inside 20' from the cluster center. Since 
a velocity was not determined for all possible candidates, some artifacts (mostly from the fibre 
allocation process) are visible in the overall distribution. 
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Fig. 4. — Velocity histograms for stars within 57' of the center of u Cen. The black line outlines 
the velocity distribution for stars with 20' < r' < 30' and the red line outlines the distribution for 
stars with 30' < r' < 57'. The systemic velocity of the cluster is 232 km s _1 and the dotted vertical 
lines indicate the velocity interval from which candidate cluster members are selected. 
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which nevertheless remains detectable above the high velocity tail of the field star distribution. 
Using the same selection interval as for the inner region, the 46 candidate cluster members have a 
mean velocity of 227 km s~ x and a velocity dispersion of 17 km s _1 . The increased velocity disper- 
sion over that seen for the inner region presumably suggests that this outer sample is contaminated 
to some extent with non-cluster members^] These stars will be referred to as the outer cluster 
sample. 

Two other criteria can be used to investigate the question of cluster membership. First, noting 
that the expected proper motion of cluster members is small (cf. Dinescu et al.|[l999a ), we can use 



published proper motion catalogs to weed out from the sample any stars with significant motions, as 
they are unlikely to be cluster members (though equally a small proper motion does not guarantee 
cluster membership). Second, while it is well established that to Cen has a substantial internal 
abundance range, we can still use the spectra to define relations between various line strength 
parameters and again seek outliers from the relations defined by the majority of stars. 

Considering first proper motions, since the majority of the candidates are too faint to be 



included in the UCAC2 catalog (Zacharias et al. 2004 1, we extracted from the SuperCOSMOS 



Science Archive^] proper motions for all the stars in the outer cluster sample (due to crowding 



effects we did not consider the proper motions for the inner cluster sample). Hambly et al. (2001 1 
indicate that for the magnitude range of our program stars, the estimated proper motion accuracy 
is of order 10 mas/yr. The mean motions in RA and Dec for the 46 stars were calculated along 
with the standard deviations, which were 6.3 mas/yr in RA and 5.5 mas/yr in dec, respectively. 
One star, 6J3J.259 (cf. Table [3]) lay outside the error ellipse defined by Sa^^a and Sa^^s- It was not 
considered in the analysis of spectral line strengths. 

For line strength parameters we first removed the overall continuum shape of the spectra by 
fitting a lower order polynomial using the IRAF routine continuum. Following, for example, |Beers 



et al. (19991, we then measured via numerical integration the equivalent widths of the Call K line 
and the G-band of CH and, via gaussian fits, the equivalent widths of the and H7 hydrogen 
lines (cf. Fig. [2]). The continuum and feature band passes used are listed in Table [lj The K- 
line and G-band strengths serve as metal abundance indicators while the hydrogen line strengths 
serve as an indicator of effective temperature. The mean uncertainties in the line strengths were 
then determined from the standard deviations of the measurements for the seven stars in the inner 
sample that have at least 6 observations. For stars with multiple spectra the individual line strength 
measurements were averaged. Only for one star was it not possible to measure all 4 indices: for 
star 8_6_15256 the H<5 line was significantly affected by noise and no reliable measurement could be 
made. The other line strengths, however, are consistent with cluster membership. 



1 We have been allocated additional time on the Anglo- Australian Telescope to measure more precise velocities for 
stars in this sample, in order to improve the constraints on the velocity dispersion in the outer parts of the cluster. 

2 http:/ /surveys. roe. ac.uk/ssa/index. html 



- 12 - 



Given that the inner sample is likely to consist largely of cluster members, we used this sample 
to define relations between the K-line strength and the average hydrogen line strength Wjj = 
[W(H)s + W(Hy)]/2, and between the G-band strength and Wh for probable cluster stars. These 
relations are shown in the upper panels of Fig. [5} Turning first to the upper left panel, we can 
see that, as expected given the well established abundance range in oj Cen, there is a significant 
range in the K-line strengths for the inner sample stars that is anti-correlated with hydrogen line 
strength. The solid line in the panel represents a least squares fit to the data, while the dotted 
lines represent ±2x the rms about the fit after excluding the four stars plotted as open circles. 
Two of these are clearly deviant, but exclusion of the two stars with the weakest K- lines deserves 



some comment. The abundance distribution in to Cen is well established (e.g. Norris et al. 19961 
and is characterized by a steep initial rise to the mean abundance followed by a long tail to higher 
abundances. Excluding the two stars with the lowest K-line strengths, the distribution of K-line 
strengths in the panel reflects this general shape and for that reason the two stars have been 
excluded. A small number of stars also lie to the right of the +2<r line in the panel. These are 
nevertheless considered members because they all show very strong G-bands and other features 
indicating probable carbon over-abundances (see Fig. [6]). The exclusion of 4 stars from the inner 
sample as possible non-members is consistent with the number of field stars expected in the adopted 
cluster velocity range, given the number of stars in the velocity intervals adjacent to that for the 
cluster seen in the solid histogram in Fig. [4| 

The upper right panel of Fig.[5]shows the distribution of G-band strengths for the inner sample. 
Once again the solid line is a least squares fit and the dotted lines are ±2x the rms about the fit 
after again excluding the four stars plotted as open circles. There is clearly a very large range in 
G-band strengths for the stars in this sample. Many of the stars with strong G-bands also show 
additional CH absorption in the region around A4370A; some examples are shown in Fig. [6j The 
existence of stars in u Cen with a range of overabundances of carbon is well established; see, for 



example, Stanford et al. (2007) and the references therein. 



The lower left panel shows the W^,Wj^ for the stars in the outer sample. The solid and dotted 
lines are reproduced from the equivalent panel for the inner sample. Given that there is evidence 
for a radial abundance gradient in the cluster, in the sense that metal-rich objects are less frequent 
in the outer regions (e.g. |Norris et aL|l996| ) , it seems unlikely that the outer sample would contain 



comparatively more metal-rich objects than the inner sample. Consequently, for the stars with 
strong K-lines, we have considered stars with Wk ^ 10.6 A to be non-members even if they are 
within the ±2<r boundaries. In many cases these stars have G-band strengths that fall outside in 
the ±2<7 boundary in the lower right panel, confirming their probable non-member status. All-in-all 
we classify 6 of the outer sample stars as non-members. As noted above there is seventh probable 
non-member in the outer sample based on a discrepant proper motion. Interestingly, like the inner 
region, the outer region also contains a number of CH-strong stars. Examples are shown in Fig. 
|6j Inspection of the dashed line histogram in Fig. [4] suggests that excluding 7 stars from the outer 
sample as likely non-members is perhaps a conservative number; there may well be an additional 
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Fig. 5. — Line strength measurements. The left panels show the equivalent width of the Call 
K line plotted against the mean equivalent width of the H<5 and H7 hydrogen lines. The right 
panels plot the G-band (CH) strength against hydrogen line strength. The upper panels are for the 
inner sample while the lower panels are for the outer sample. Open circles represent stars that are 
classified as probable non-members of the cluster. The solid lines are least squares fits to the inner 
sample, with the dotted lines representing ±2x the rms about the fit . The fits are reproduced in 
the lower panels. Mean error bars are shown in the bottom left corner of each panel. 
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number of approximately the same size of non-member stars in this sample that we have not been 
able to conclusively identify. Such identification would require velocities with a higher degree of 
precision than that presented here as well as knowledge of the cluster velocity dispersion in this 
outer region. 

One star deserves special mention. This is star 10_2_28 whose spectrum is shown in the bottom 
panel of Fig. [6| With strong C2 features at A ~ 4700A, strong CH absorption features in the region 
AA4250-4430A and strong CN absorption features at A ra 3883A and A ss 42 15 A, this star is 



clearly a carbon-star. Such stars are not unknown in oj Cen (e.g. van Loon et al. 20071 but the 



known members stars with similar spectra, such as ROA 577 (see Fig. 25 of van Loon et al 



20071, are found on the red giant branch at considerably brighter magnitudes. Determining the 
radial velocity of 10_2_28 is not straightforward given the very different spectral type. Nevertheless, 
cross-correlations with the latest spectral-type among the observed radial velocity standards and 
with probable oj Cen members such as 8_8_4146 (cf. Fig. [6J indicate a velocity of ~260 db 15 km 
s _1 . This is nominally within the velocity range of probable cluster members but the star lies at a 
distance of 47.3' from the cluster center. Given that the number of cluster members at this radius 
is small, and that cluster members of this type are comparatively rare, we conclude that this star 
is likely not a member of the cluster. Interestingly there is a second carbon star that lies in the 



field of the cluster but which is clearly not a cluster member (ROA 153, see Smith & Wing 19731. 



In Table[2]we list the ID number, the radial distance from the cluster center, the right ascension 
and declination, the number of observations, the radial velocity and the V and V — I photometry 



(from Coleman Sz Da Costa 20081 for the 154 probable oj Cen members categorized in this work. 
At a distance of 54.7' from the cluster center, star 9_4_1918 is the most radially distant probable 
cluster member inside the adopted tidal radius; there are 9 such stars beyond 40' from the cluster 
center. For completeness in Table [3] we list the same information for the 12 stars in the inner 
and outer cluster samples that are classified as probable non-members, including the carbon star 
10_2_28, despite having radial velocities comparable to that of the cluster. Similarly, Table [4] lists 
the same information for the 1183 stars inside the adopted cluster tidal radius that have velocities 
outside the velocity interval (190-270 km s _1 ) used to select cluster member candidates. 

Fig. [7] shows the photometry for the probable members superposed on the photometry for all 
stars in the 10'-25' radius region (cf. Fig. [I]). As expected given the cluster metallicity distribution, 
the blue side of the color distribution of the probable members is sharply defined but there is a 
considerable spread to redder colors. 



3.2. Stars beyond the tidal radius 

In Fig. [8] we show the velocity distribution for stars beyond the adopted tidal radius. The 
distribution is clearly dominated by disk stars which have heliocentric velocities near zero (for 
reference, zero velocity in the galactic rest frame at the Sun's location corresponds to a heliocentric 
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Fig. 6. — Spectra of stars showing enhanced carbon features at the G-band, A « 4300A, at addi- 
tional A-X CH bands, A « 4370A, and at C 2 bands, A « 4700A. The upper two stars come from 
the inner sample while the third is from the outer sample. The star in the bottom panel is clearly 
a carbon star. Its membership status is discussed in the text. 
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Fig. 7. — The V, J photometry for the 154 probable members of u Cen superposed on the photom- 
etry from Fig. [TJ 
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velocity of 168 km s _1 in the direction of lo Cen, assuming a galactic rotation velocity of 220 km 
s" 1 ). The insert in the figure shows a close-up of the distribution in the vicinity of the velocity 
of the cluster. While the numbers of stars are small and therefore subject to significant statistical 
fluctuation, there is no obvious indication of an excess of stars at or near the cluster velocity. 

Nevertheless there are 28 stars in the ±40 km s _1 velocity range shown in the insert frame 
in Fig. [8| Of these, two have proper motions inconsistent with association with lo Cen and all 
(including the stars with discrepant proper motions) but 9 occupy locations in the (Wk,Wh) arid 
(WcWff) diagrams that are inconsistent with lo Cen membership (cf. Fig. [5j. Most of these likely 
non-member stars have relatively large values of Wx; 18 have Wk > 9.5A (cf. Fig. [5]). They are 
most likely members of the Galaxy's bulge population lying at distances beyond that of lo Cen. 
The line-of-sight through the cluster has a minimum distance of 10.4 kpc from the Galactic Center 
at a height of 3.4 kpc above the Galactic plane. At this distance the stars would have My ~ 
+0.5, approximately the magnitude of the red clump, which is not inconsistent with the proposed 
interpretation. 

For the 9 stars with line strengths consistent with cluster membership, two have velocities of 
192 and 195 km s , respectively, while a third has V r = 261 km s . The remaining six have a 
mean velocity of 233 km s -1 , corresponding closely to that for the cluster, and a dispersion of only 9 
km s _1 , which given the velocity errors, is consistent with an intrinsic dispersion of near zero. These 
are exactly the characteristics expected for an extra-tidal population and we tentatively identify 
these six stars as likely extra-tidal (former) cluster members. Table [5] lists their properties. The 
spectrum of one of the stars is shown in the bottom panel of Fig. [2j Given the small number and 
the shape of the region sampled (cf. Fig.[3|, no inferences can be made from the spatial distribution 
of these stars, except to note that the furthest lies almost exactly twice the adopted tidal radius 
from the cluster center (which is also the approximate limit of the area surveyed) and that the 
surface density corresponds to 1.2 stars per square degree for the region between one and two tidal 
radii. Clearly more accurate velocities and abundance information are required to more firmly 
establish the association or otherwise of these stars with lo Cen, but for the present we will assume 
that they are all extra-tidal members of the cluster and treat derived quantities as upper limits. 
The alternative assumption, that all or some of these stars are simply members of the field halo 
population, cannot be ruled out. 

The immediate question to ask is what fraction of the mass of lo Cen is potentially represented 
by the extra-tidal stars. Assuming that the extra-tidal population is identical with that of the 
cluster, we can estimate the fraction as follows. Using a core radius of 2.4' and a tidal radius of 57' 



and the integral formulae given in King ( 1962 1, the percentage of the total luminosity of the cluster 
contained between 20' and the tidal radius is approximately 10 per cent. Assuming mass follows 
light, this figure also represents the fraction of the total mass contained in this radial region. The 
region is also characterized by the 154 likely members we have identified, whereas using the surface 
density calculated above, the region between one and two tidal radii is estimated to contain, at most, 
approximately 10 extra-tidal members with the same characteristics. Consequently, combining 



- 18 - 




-100 100 200 300 400 



Heliocentric Velocity 

Fig. 8. — Velocity histogram for stars lying further than 57', the adopted tidal radius, from the 
center of u Cen. The insert panel shows the distribution near the heliocentric velocity of the cluster 
and the dotted lines outline a ±40km s^ 1 window about that value (cf. Fig. EJ. 
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these numbers, we conclude that an upper limit on the amount of mass in extra-tidal material lying 
between one and two tidal radii from the cluster centre is of order 0.7 ± 0.2 per cent of the total 
cluster mass. It may well be considerably less. The uncertainty quoted comes solely from assuming 
Poisson statistics apply to the sample of 6 identified extra-tidal member candidates. This limit is 
comparable to the results of Kiss et al.| (2007) who estimated a ~1% upper limit on the extra-tidal 



populations of three globular clusters (NGC 288, M30 and M55) but differs from that of Grillmair 



& Johnson (20061 who estimate that at least 3% of the present cluster mass is associated with the 



tidal tails of the globular cluster NGC 5466. The biggest contrast is of course with Pal 5, where 



the mass in the extensive tidal tails exceeds the mass of the cluster itself (e.g. Odenkirchen et al. 
20031. 



3.3. The surface density profile 



In Fig.[9]we show the (circularly averaged) surface density profile of u Cen from the center to the 
outermost regions, covering a factor of ~10 5 in surface density. The profile combines the surface 
photometry and star-count data from Da Costa (1979)]^ together with the circularly averaged 



surface density data from Coleman & Da Costa (2008) and that derived from our cluster member 



sample. The points from Coleman & Da Costa (2008) are for stars from just below the turnoff and 



brighter (cf. Fig. [TJ and are for regions where image crowding is not a concern. The background 
density was determined from the extensive region beyond the tidal radius covered in the imaging 
survey. The stars included correspond to those that dominate the integrated light and thus no 
systematic effects are expected in scaling the surface density count data to the surface brightness 
observations (cf. |Da Costa 1979) . In any case the long 2-body relaxation time of u Cen (e.g. 



Harris 1996), especially in the outer parts, means that significant mass segregation in this cluster 



is unlikely. Indeed van de Ven et al. (2006) find no dynamical evidence for any change in M/L with 
radius, as expected given the relatively long relaxation time. 

The surface density points from the cluster member sample were determined as follows. First, 
we grouped the radial distribution of the 154 probable oj Cen members into annuli. Mean surface 
densities were then calculated with errors determined from the number of stars in each annulus. 
While Fig. [3] shows that each annulus has not been completely surveyed, we assume that there 
is no radial variation in the degree of incompleteness. The resulting surface densities were then 
matched with the existing data in the region between 20' and 26'. As is evident from Fig. [9] the 
agreement between the new and existing data at these, and larger radii, is excellent. This serves 
as confirmation that our sample, while not complete, is suitably representative. Overall the profile 
shows a rapidly declining surface density distribution with a limiting radius consistent with the 
value we have adopted for the cluster tidal radius, shown by the vertical line in the Figure. Fig. 



3 These data differ insignificantly from those of Meylan ( 1987 1 who adopted a very slightly different normalization 
of the star count data to the surface photometry. 
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[9] also shows the density inferred for the extra-tidal population in the region between one and two 
tidal radii, if the six candidates identified here are all associated with the cluster. 

Given the extensive range in surface density, the complete profile shown in Fig. [9j not sur- 
prisingly, cannot be adequately represented by analytic formulae such as those of King| ( |1962 1 , or 



by simple models (e.g. King 1966). However, more complex models, for example those of Meylan 



(1987), Meylan et al. (1995) and Merritt et al. (1997), do appear capable of fitting the entire profile 



We note in passing that the models of van de Ven et al. (2006) adopt a tidal radius of 45', which 



on the basis of Fig. [9j is somewhat too small. The core radius of the cluster is well established at 
r c « 2.5' (e.g. |King|1962| |Da Costa||1979| fTrager et al.|[l995| ). Table © lists the surface density data 



from Coleman & Da Costa (2008), and from this work, in a manner that is directly compatible 



with Table 1 of Meylan (1987). 



Discussion 



The photometric work of Law et al. ( 2003| ) (see also Coleman & Da Costa 2008) concluded 
that u Cen lacks any notable extra-tidal structure, with the extensive tidal tails originally proposed 



by Leon et al. (2000) arising from the lack of allowance for the effects of reddening variations across 



the field (a possibility that Leon et al. (2000) themselves acknowledged). Our spectroscopy based 



results are consistent with, and strengthen, this result. 

As noted in §1, the presence of extra-tidal structure, particularly in the form of tidal tails, is 



relatively well established in a number of globular clusters (e.g. Grillmair et al. 1995 Leon et al. 



2000 Grillmair & Johnson 2006 and the references therein) with the most striking example being 



Pal 5 (Odenkirchen et al. 2003 Grillmair & Dionatos 2006). Clearly, however, oj Cen does not 



obviously fall into this category. The most likely explanation for this is as follows. The largest 
dynamical influence on co Cen is the 'tidal shocking' that occurs each time the cluster passes through 



the Galactic plane (e.g. van de Ven et al. 2006), as is also the case for Pal 5 (cf. Dehnen et al. 



2004). van de Ven et al. (2006) present a detailed model of the structure of u> Cen, constrained 



by high quality proper motions and radial velocities for many thousands of cluster member stars, 
and we can use this model to estimate the effects of tidal shocks on the outer parts of the cluster. 



With the parameters given in §9.5 of van de Ven et al. (2006), which rely on the 10 Cen orbital 



parameters listed by Dinescu et al. (1999b), the average change in the velocity of an u Cen star as 



a result of the impulse that occurs as the cluster traverses the Galactic plane is estimated as \Av\ 
~ 0.17 r' km s , where r' is the radial distance of the star from the cluster center in arcmin. This 
is to be compared with the local escape velocity, which can be approximated by twice the local 



velocity dispersion (e.g. Binney & Tremaine 1987 p. 490). The form of the velocity dispersion in 



the outer parts of ui Cen is not well established (even the observations of Scarpa et al. (2003) reach 
only a maximum of ~30' from the center) but it is clear that it is only in the outermost parts of 
the cluster that the tidal shock effects are large enough to add sufficient energy to (some) cluster 



stars to cause them to escape. Indeed as van de Ven et al. (2006) note, the phase-space structure of 
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Fig. 9. — The (circularly averaged) surface density profile for u Cen. Open squares and star 
symbols are surface photometry and star-count data, respectively, from Da Costa (19791. Filled 
triangles are data from Coleman & Da Costa (2008) while the filled circles are derived from the 
cluster member sample determined in this paper. The vertical scale corresponds to the surface 
density of stars in the cluster member sample. It can be converted to the surface brightness scale 
of Da Costa (1979) by subtracting 8.318. The vertical dot-dash line indicates the adopted tidal 
radius. The horizontal dashed line shows the surface density inferred for the region between 1 and 
2 tidal radii if the six candidate extra-tidal stars are indeed associated with the cluster. The dotted 
lines represent the statistical uncertainty in that density. 
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the outer parts of oj Cen may well be completely dominated by the effects of the relatively frequent 



tidal shocks: van de Ven et al. (20061 estimate that the cluster is immersed in the disk and feels 
the additional gravitational field for ~10% of the ~120 Myr orbital period. The main point here 
though is that the region of oj Cen which is most strongly influenced by the tidal shocks contains 



only a very small fraction of the total mass of the cluster, unlike the situation for Pal 5 (Dehnen 
et al. 2004). For example, only ~1% of the oj Cen mass lies between ~40' and the tidal radius. 



Consequently, it is likely that tidal shocks are indeed driving the escape of stars consistent with 
the limits on the extra-tidal population observed here. Further, we note that if stars leave the 
cluster with a relative velocity of ~1 km s , then that is sufficient for the star to traverse the 
region between 1 and 2 tidal radii in less than the orbital period of the cluster around the Galaxy. 
Thus even if tidal shocks are driving the escape of a small number of stars, the number remaining 
close to the cluster will be small, again consistent with the observed lack of extra-tidal stars in our 
study. Indeed this suggests that if it were possible to isolate putative oj Cen main sequence stars at 
and beyond the tidal radius against the enormous field contamination, it may well be possible to 
characterize tidal tails from oj Cen. Such discrimination might be possible with GAIA data though 
many square degrees of sky would need to be searched. 

In the context of the model in which oj Cen is the remnant nucleus of a disrupted former 
dwarf galaxy, it is clear from the above that any stars which are more loosely bound than those 
in the outer parts of the remnant nucleus/star cluster will be even more strongly affected by the 
tidal shocks. Consequently, they would have been lost from the vicinity of the remnant nucleus/star 
cluster at relatively early stages in the interaction, once the disrupting dwarf galaxy has approached 



the current orbit of oj Cen. The dynamical models of Bekki Sz Freeman (2003) and Tsuchiya et al. 



(2003 ) show that this occurs within about 3 Gyr of the initiation of the interaction. Given that any 



gas in the original dwarf would have been used up in star formation (cf. Bekki & Freeman 2003 ) 
(and any remaining stripped out) in the first phases of the interaction, and that oj Cen contains 
only stars with ages in excess of ~10 Gyr, (e.g. |Sollima et al.|2005[ |Stanford et al.|2006| ) it is evident 
that the more loosely bound stars have been lost from the vicinity of oj Cen many Gyr ago, and 
are now widely dispersed around the Galaxy. The kinematic signatures of local field stars that may 
have had their origin with the disrupted system that once contained oj Cen have been modeled, 
and such stars may well have been observationally detected (e.g. Dinescu [2002 Mizutani, Chiba & 



Sakamoto 2003; Chiba & Mizutani||2004 Meza et al. 2005) 



In this sense the currently disrupting dwarf galaxy Sagittarius and its central massive star 
cluster M54 provide a similar example to the postulated scenario for oj Cen, but on a longer 
timescale - eventually the Sgr field stars will disperse uniformly around the galaxy but like oj Cen, 
M54 will retain its dynamical integrity for a very long time. 
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5. Conclusions 

We have carried out an extensive radial velocity and line-strength analysis of observations of 
a large sample of stars that covers the outer parts of the globular cluster uj Cen and a substantial 
area beyond the cluster's tidal radius. While the analysis has led to the establishment of cluster 
membership for ~150 stars in the outer parts of the cluster, and revealed that the steeply declining 
radial surface density distribution of these stars is consistent with other extant data, the survey 
has identified only a very small number of objects that can be plausibly classified as extra-tidal 
cluster stars. The upper limit on the amount of mass potentially contained in the region between 
1 and 2 tidal radii is less than one percent of the present cluster mass. We conclude therefore that 
there is no compelling evidence for any significant extra-tidal population associated with uj Cen. 

The authors are grateful for the excellent support at the telescope provided by the staff of the 
Anglo- Australian Observatory. The comments of the anonymous referee led to improvements in 
the clarity of the paper. This research has made use of the SIMBAD database, operated at CDS, 
Strasbourg, France. 

Facilities: AAT (2dF), ANU Siding Spring Observatory lm telescope 
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Table 1. Spectral feature and continuum bandpasses 



Feature 


Feature Bandpass 


Blue Continuum 


Red Continuum 




(A) 


Bandpass (A) 


Bandpass (A) 


Call K 


3925.-3943. 


3909.-3921. 


4010.-4030. 


m 


4093.-4109. 


4010.-4030. 


4120.-4131. 


G-band (CH) 


4294.-4316. 


4230.-4248. 


4415.-4435. 


H 7 


4334.-4347. 


4230.-4248. 


4415.-4435. 



Table 2. Probable cluster member stars 



ID 


R 


Velocity 


N 


RA 


Dec 


V 


V -I 




(arcmin) 


(km s _1 ) 




(J2000) 


(J2000) 






8_4_3206 


20.0 


239 


1 


13 27 47.3 


-47 45 44 


16.42 


0.97 


8_7_15831 


20.1 


224 


1 


13 25 06.5 


-47 17 42 


15.80 


0.96 


8_6_16385 


20.1 


219 


1 


13 24 56.3 


-47 36 35 


16.39 


0.93 


8_8_3134 


20.2 


221 


2 


13 25 47.7 


-47 11 00 


16.63 


0.89 


8_6_15256 


20.2 


197 


1 


13 25 11.4 


-47 41 03 


16.66 


0.90 


8_8_2219 


20.2 


228 


1 


13 26 04.4 


-47 09 40 


16.42 


1.02 


8_4_4900 


20.3 


248 


1 


13 27 20.0 


-47 48 03 


16.23 


0.97 


8_2_1336 


20.3 


253 


1 


13 28 43.3 


-47 24 38 


16.19 


1.00 


7_3_261 


20.3 


222 


1 


13 24 46.0 


-47 30 14 


16.68 


1.07 


5_4_1514 


20.4 


209 


1 


13 27 09.5 


-47 08 34 


16.47 


0.93 



Note. — Table [2] is published in its entirety in the electronic addition of the 
Astronomical Journal. A portion is shown here for guidance regarding its form 
and content. 
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Table 3. Probable non-cluster-member stars 



ID 


R Velocity N RA 


Dec 


V V-I 




(arcmin) (km s" 1 ) (J2000) 


(J2000) 





7_2_312 


20.7 


197 


2 


13 


24 


44.4 


-47 


25 


56 


16.21 


0. 


.96 


11_1_4329 


21.4 


226 


1 


13 


27 


27.7 


-47 


48 


47 


16.53 


0. 


.96 


8_1_1381 


22.5 


232 


1 


13 


28 


24.3 


-47 


13 


33 


16.55 


0. 


.88 


11_8_1684 


23.0 


232 


1 


13 


26 


00.2 


-47 


50 


14 


16.06 


1 


.04 


11_1_3775 


32.5 


219 


1 


13 


27 


39.0 


-47 


59 


55 


15.56 


1 


.22 


5_6_92 


34.9 


202 


1 


13 


26 


36.1 


-46 


53 


48 


16.06 


1. 


.03 


12_8_1988 


37.9 


206 


1 


13 


29 


30.5 


-47 


54 


25 


16.35 


1 


.15 


6_6_1259 


45.0 


197 


2 


13 


29 


56.3 


-46 


57 


16 


16.73 


0. 


.89 


10_2_28 


47.3 


260 


1 


13 


24 


45.8 


-48 


11 


24 


16.17 


1 


.14 


10_1_2777 


48.0 


254 


1 


13 


22 


31.9 


-47 


50 


17 


16.43 


0. 


.88 


4_3_2022 


54.6 


256 


2 


13 


23 


07.1 


-46 


48 


42 


16.15 


1 


.21 


11_3_191 


55.1 


243 


1 


13 


29 


17.7 


-48 


17 


31 


15.54 


1 


.16 



Note. - Star 6_6_1259 excluded on basis of discrepant proper motion; all 
others excluded on the basis of discrepant line strengths; star 10_2_28 is a carbon 
star. 
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Table 4. Non-member stars 



ID 


R 


Velocity 


N 


RA 


Dec 


V 


V — I 




(arcmin) 


(km s _1 ) 




(J2000) 


(J2000) 






11_8_135 


20.0 


-11 


2 


13 26 44.3 


-47 48 38 


16.05 


1.02 


8_4_4179 


20.2 


-42 


1 


13 27 31.3 


-47 47 16 


16.43 


1.02 


11_1_5863 


20.2 


34 


10 


13 26 59.0 


-47 48 44 


16.23 


1.08 


8_8_1220 


20.5 


-5 


8 


13 26 21.6 


-47 08 30 


16.39 


0.85 


8_1_1736 


20.6 


-11 


1 


13 28 14.9 


-47 14 39 


16.24 


0.94 


8_2_2607 


20.6 


168 


1 


13 28 23.8 


-47 16 18 


15.67 


1.05 


8_8_2721 


20.6 


26 


1 


13 25 55.6 


-47 09 52 


16.65 


0.95 


8_5_4726 


20.7 


6 


1 


13 25 17.3 


-47 43 02 


16.29 


1.03 


8_8_4580 


20.7 


-19 


1 


13 25 12.7 


-47 15 09 


16.64 


0.97 


5_4_1370 


20.9 


55 


8 


13 27 19.0 


-47 08 30 


16.19 


0.99 



Note. — Table [4] is published in its entirety in the electronic addition of the 
Astronomical Journal. A portion is shown here for guidance regarding its form 
and content. 



Table 5. Candidate extra-tidal cluster members 



ID 


R 


Velocity 


N 


RA 


Dec 


V 


V -I 




(arcmin) 


(km s _1 ) 




(J2000) 


(J2000) 






4_7_288 


66.5 


238 


1 


13 22 03.5 


-46 42 37 


16.43 


0.98 


7_8_1878 


70.1 


238 


1 


13 20 07.1 


-47 09 50 


15.70 


1.28 


14_3_2302 


93.0 


246 


1 


13 27 49.7 


-49 01 00 


16.57 


1.08 


1_2_745 


99.0 


230 


1 


13 23 42.4 


-45 54 46 


15.90 


1.18 


15_5_941 


112.6 


219 


1 


13 30 36.3 


-49 14 28 


16.50 


1.02 


3_8_141 


113.3 


227 


1 


13 31 01.3 


-45 44 08 


16.42 


1.01 
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Table 6. u Cen surface density data 



log r 


log 


SD 


log r 


log SD 


[arcmin] 


[10.0 V mag/arcmin 2 ] 


[arcmin] 


[10.0 V mag/arcmin 2 ] 


1.115 


-2.076 


± 0.011 


1.323 


-2.769 ± 0.063 


1.177 


-2.308 


± 0.014 


1.362 


-3.022 ± 0.080 


1.231 


-2.449 


± 0.016 


1.398 


-3.165 ± 0.090 


1.279 


-2.674 


± 0.022 


1.448 


-3.806 ± 0.129 


1.323 


-2.809 


± 0.026 


1.513 


-3.871 ± 0.115 


1.362 


-2.928 


± 0.031 


1.575 


-3.933 ± 0.115 


1.398 


-3.090 


± 0.040 


1.629 


-4.561± 0.239 


1.432 


-3.164 


± 0.044 


1.677 


-4.734 ± 0.286 


1.463 


-3.343 


± 0.061 


1.729 


-5.108 ± 0.383 


1.492 


-3.422 


± 0.069 






1.532 


-3.766 


± 0.097 






1.589 


-4.019 


± 0.160 







Note. — The first 2 columns are from Coleman & Da Costa (20081, 
the remainder from the cluster member sample determined here. The 



surface densities have the same normalization as Table 1 of Meylan 



(1987|. Add 0.832 (equivalent to a central surface brightness of 16.81 
V mag/arcsec 2 ) for log(surface density) values in the adopted units. 



